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Abstract 

Consider a closed lipid membrane (vesicle), modeled as a two-dimensional surface, described by a 
geometrical hamiltonian that depends on its extrinsic curvature. The vanishing of its first variation 
determines the equilibrium configurations for the system. In this paper, we examine the second varia- 
tion of the hamiltonian about any given equilibrium, using an explicitly surface covariant geometrical 
approach. We identify the operator which determines the stability of equilibrium configurations. 

1 Introduction 

Many of the physical properties of a lipid membrane (or vesicle) are captured by a hamiltonian which 
describes the degrees of freedom of an idealized two-dimensional membrane surface. Typically, this hamil- 
tonian is a sum of terms each of which either penalizes or constrains some geometrical characteristic of this 
surface. A term quadratic in the mean extrinsic curvature is a measure of the energy penalty associated 
with the bending of the membrane [Q ^, [|. This term is supplemented by various constraints: the area 
as well as the enclosed volume are usually fixed. A constraint or (a term in the energy) linear in the mean 
extrinsic curvature describes an asymmetry between the lipid bilayers which constitute the membrane [Q . 
(For a review see |, |, |, |, g.) 

In equilibrium, this hamiltonian will be stationary with respect to arbitrary infinitesimal deformations 
in the surface. This infinitesimal deformation has both a tangential part and a normal part. The former 
corresponds to a reparametrization of the surface, and so can contribute only in boundary terms. The 
remaining normal infinitesimal deformation will completely describe the physical state of deformation 
of the membrane in its bulk. As a consequence, the equilibrium condition is a single equation, known 
somewhat prosaically as 'the shape equation' JTofl , In the thirty or so years since the introduction of the 
model, a lot has been learnt about the solutions of this equation, both analytically and numerically (see, 
for example, Q and the references it contains). 

In this paper, we will examine the second variation of the hamiltonian. From one point of view, the 
calculation of the second variation about an equilibrium configuration is necessary in order to assess its 
degree of stability; in addition, the fluctuations of the membrane due to its interaction with its environment, 
are encoded completely at lowest order in the second variation of the hamiltonian. This problem has, of 
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course, been addressed before and in both contexts. In the former, for example, perturbations about 
spherical or cylindrical configurations have been examined analytically jn], [l3| |TJ] . A decomposition of 
geometries with a convex hull into spherical harmonics has also proven very effective [fl5f . In the statistical 
mechanical context, while arbitrary background geometries have been considered, the focus has been on 
the part of the second variation contributing at one-loop to the renormalization of the parameters of the 
theory || 0, |J, In any case, there is a gap in perturbation theory which we would like to fill. 
We will write down the second variation of the hamiltonian about an arbitrary equilibrium configuration, 
approaching the problem from a manifestly geometrically covariant point of view. Two features of the 
approach we will adopt should be emphasized. At the level of the second variation, it is no longer justified 
to ignore tangential deformations of the membrane. This is because a finite tangential deformation unlike 
its infinitesimal counterpart is not a simple reparametrization of the surface. Indeed, it was by examining 
finite tangential deformations of a given surface of constant negative Gaussian curvature preserving the 
curvature that Backlund, Bianchi and others generated new surfaces pp[ ]. We show, however, that when 
the background geometry is an equilibrium, tangential deformations contribute only boundary terms to 
the second variation of the Hamiltonian, and therefore for the consideration of bulk fluctuations about 
equilibrium, they can be neglected without incurring any error. A second simple but subtle technical 
point is the exploitation of the fact that the variation of the divergence of a vector density is equal to 
the divergence of the variation. Though this is not so essential computationally at the quadratic level we 
will work at, by facilitating the isolation of boundary terms, it does make higher order expansions in the 
fluctuations about a non-trivial background feasible in practice. 

The paper is organized as follows. In section 2, we introduce our conventions for the geometrical 
description of a hypersurface embedded in R N+1 . We extend our considerations to an A-dimensional 
hypersurface because of the little extra cost involved and because of potential applications in other geo- 
metrical problems in soft matter physics. In section 3, we examine how the geometry changes under an 
infinitesimal deformation of the hypersurface. This is used in section 3 to derive the first variation of a 
geometrical hamiltonian for lipid membranes, and to obtain its Euler-Lagrange derivative. In section 4 
the the second variation is examined. In particular, we check that it vanishes at equilibrium when the 
deformation corresponds to a rigid normal translation. We end in section 6 with a few concluding remarks. 



2 Geometry 

We begin by describing briefly the geometry of a hypersurface embedded in Euclidean space R N+1 . This 
allows us to introduce our conventions. We will emphasize the peculiarities associated with a surface 
embedded in R 3 . At the end of this section, we also address the issue of identifying the low order inde- 
pendent reparametrization invariants one can construct from the geometrical quantities that characterize 
the hypersurface. 

Consider an orientable hypersurface S embedded in R N+1 . This surface can be specified locally in 
parametric form by N + 1 shape functions, 

x = x(n, (i) 

where x ~ x^ L = (x 1 , ■ ■ ■ , x N+1 ) are coordinates for R N+1 , £ a arbitrary coordinates on the surface £ 
(a, b, ■ ■ ■ = 1, • • • , A), and X = (A 1 , • • • , X N+1 ) are the shape functions. 
The Euclidean metric on R N+1 induces the metric g a b on S defined by, 

g a b ■= e a ■ e b , (2) 

where the N tangent vectors are defined by e a (£ a ) = <9 a X (d a := d/d^ a ). Latin indices are lowered and 
raised with g ao , and its inverse g ab , respectively. The metric g ao determines the intrinsic geometry of the 
hypersurface X. It defines the unique torsionless covariant derivative V a compatible with it, i.e. satisfying 
V a <7& c = and (V a Vb — Vf,V a )/(^ a ) = for some surface function /(£"). In terms of the Christoffel 
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symbol r^ b , , acting on a hypersurface vector V a , it reads 

V a V b = 8 a V b + T b ac V c , (3) 

where ^ 

Kb ■= 9 cd e d ■ d a e b = -g cd (d a g bd + d b g ad - d d g ab ) . (4) 

The intrinsic Riemann curvature lZ a bcd of V a is defined by 

(V a V fc - V fc V a )V c =: TZ c dab V d . (5) 

In terms of the Christoffel symbol, the Riemann curvature tensor is given by 

Tl a bc d = d c T^ b - d d T a cb + T a ce T e db - T^ e T e cb . (6) 

The Ricci tensor is defined by contraction, lZ ab ■= 7^ c ac &; the scalar curvature 71 is defined by 1Z :— g ab lZ ab . 

When the hypersurface S is two-dimensional, N = 2, the Riemann curvature tensor is completely 
determined by the scalar curvature: 

-JZ 

TZabd = -^{gac9b d - 9a d 9bc) ■ (7) 

Note that, as a consequence, we then have that in two dimensions the Einstein tensor vanishes, 

Gab ■= Tlab - \ll9ab = . (8) 

The scalar curvature of a two-dimensional surface is related to the Gaussian curvature G of the surface by 
K = 2G. 

The simplest geometrical quantity invariant under reparameterizations of the surface one can construct 
out of the intrinsic geometry of the hypersurface X is its area, 

A:= J dA = J (9) 

where g denotes the determinant of the metric g ab . The next order invariant depending only on the intrinsic 
geometry of the surface is the average scalar curvature over the surface, J dA 1Z. For a two-dimensional 
surface with no boundary, by the well-known Gauss-Bonnet theorem, this is not only a reparameterization 
invariant, it is also a topological invariant, with 



dA TZ = 4tt(1 - g) , (10) 

where g is the genus of the surface. 

Let us turn now to the extrinsic geometry of S. The single normal vector n(£ a ) to S in R N can be 
defined in implicit form by, 

e a -n = 0, (11) 

with the normalization, 

n-n=l. (12) 

Note that these equations determine n only up to a sign: the normal can be inward or outward. Since 
the surface is assumed to be orientable, we can pick one sign consistently, and we choose the normal to be 
outward. The unit normal vector field can also be given explicitly as 

1 oi— ajv-i ( p Pi )...( P PN-i \ (io\ 



Nl^/g 
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where s^...^ and £ Q i - a N-i are the totally antisymmetric Levi-Civita symbols for R N+1 and E, respec- 
tively (ei-.-n = +1)- The factor of y/g is necessary in order to make n a scalar under reparametrizations. 
When E is closed, as we are assuming in this paper, the total volume enclosed by the surface E in 

R N+1 

is an invariant. The normal to E, as given by (|l3|), allows us to provide an alternative definition of 
the volume occupied by the interior of E as the surface integral 

V = — !— / dAii-X. (14) 

The space vectors {n, e a } form a basis adapted to the hypersurface. Their gradients along the surface 
arc themselves space vectors, and can be decomposed in turn with respect to this basis. These decompo- 
sitions constitute the classical Gauss- Weingarten equations, 

d a e b = T c ab e c - K ab n , (15) 
d a n = K ab g bc e c . (16) 

Here T c ab is the Christoffel symbol defined in (||) . The extrinsic curvature of E is given by the symmetric 
rank two surface tensor, 

K ab := n ■ d a e b = K ba . (17) 

(Note that many authors differ by a sign in this definition.) We define its trace with respect to the intrinsic 
metric, 

K :=g ab K ab =:2H, (18) 

where H represents the mean extrinsic curvature of the surface. 

The intrinsic and the extrinsic geometries of E, determined respectively by g ab and by K ab , cannot 
be specified independently. They are related by the well-known integrability conditions of Gauss-Codazzi, 
and Codazzi-Mainardi, given respectively by, 

K abcd - K ac K bd + K ad K bc = 0, (19) 
V a K bc - V b K ac = 0. (20) 

These equations follow as integrability conditions from taking a gradient along the surface of of the Gauss- 
Weingarten equations ([l5]), (0), and then the appropriate anti-symmetric part. 

The fundamental theorem for surfaces states that, given g ab and K ab , these equations are not only 
necessary, but also sufficient for the existence of an embedding with these quantities as intrinsic metric 
and extrinsic curvature. Furthermore, the embedding is unique, up to rigid motions in the ambient space 



(see e.g. |21 ) 



Contraction of the Gauss-Codazzi-Mainardi equations with the contravariant intrinsic metric g results 

K ab -KK ab + K ac K b c = 0, (21) 

K-K 2 + K ab K ab = 0, (22) 

V a K b a -V b K = 0. (23) 



For a two-dimensional surface, the contracted equations (|22| ) and ( [2 3D possess the same content as the full 
Gauss-Codazzi-Mainardi equations, ( P~9| ) and (pfj|). 

The lowest order geometrical invariant which involves the extrinsic geometry of the surface is the trace 
of the extrinsic curvature integrated over of the surface, 

M := / dAK. (24) 
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This invariant, linear in n, depends on the orientation of the surface. Note that [K] = L 1 , so that 
[M] = L. The next order reparameterization invariants are quadratic in the extrinsic curvature (w L° ), 

dAK 2 , J dA K ab K ab , 

together with the Gaussian term, J dA 1Z. In general, these three invariants are not independent, as 
follows from the contracted Gauss equation, . Furthermore, in two dimensions, the Gaussian term is 
the Gauss-Bonnet topological invariant, see (|10|). For a two-dimensional surface, at this order there is then 
only one independent local scalar associated with the embedding, and it is customary to choose J dA K 2 . 
The gradients of the extrinsic curvature V a K ab which appear in the Codazzi-Mainardi equations do not 
feature. Differential bending can be ignored up to this order. 

Let us now construct explicitly the geometrically independent terms of order L~ 2 and lower. In general, 
at order L~~ 4 , we have three independent scalars. Of the set, K 3 , KK ab K ab , K ab K bc K c a , KU, and K ab Q ab , 
the Gauss-Codazzi equations can be exploited to produce three that are independent 



J dAK 3 , JdAKK, and JdAK ab g ab . 



The last of the three vanishes identically for two dimensional surfaces. 

At order L~ 2 using the Gauss-Codazzi equation ([l9]), it is easy to show that the independent scalars 
are TZ abcd TZ abcd , TZ ab TZ ab , II 2 , K 4 , TZK 2 , and G a bK ab . Of these, for a two-dimensional surface, only II 2 , 
K 4 , UK 2 survive. Moreover, at this order gradients of the extrinsic geometry enter. We can reduce all 
terms quadratic in V a K bc to the form V a KV a K plus terms that we have already considered. We exploit 
the uncontracted Codazzi-Mainardi integrability condition ( ^p| ) to write 

J dA {V a K bc )(V a K bc ) = J dA {V a K bc ){V b K ac ) = - j ' dA K bc V a V b K ac , 

where we have integrated by parts. We now note that 

[V a , V„]K ac = K a dab K dc + TZ c dab K ad , 

so that 

J dA {V a K bc ){V a K bc ) = -J dAK b {V b V a K ac + TZ db K dc + TZ c dab K ad ) 



J dA (K bc X7 b X7 c K + TZ db K dc K b c + K cdab K ad K bc ) , 



where we have used the contracted Codazzi-Mainardi equation (|23|) in the first term. We now integrate 
by parts again the first term, and use the Gauss-Codazzi-Mainardi equations to obtain 

J dA {V a K bc ){V a K bc ) = J dA(V a K){V a K)+K ab K ab - KK ab K ab + ^TZ abcd TZ abcd . 

We conclude that for a two-dimensional surface, at this order there are four independent scalars, 

J dA K 2 , JdAKK 2 , J dA K 4 , and J dA (V a K){V a K) . 

These higher order terms appear in geometric models for the so called egg-carton membranes [^2| , and for 
tubular structures |2^]. Note that if K is treated as a simple scalar field, <f>, then the model described at 
this order is a A</> 4 theory non-minimally coupled to the intrinsic curvature. The Euler-Lagrangc equations 
of the two models are, of course, different. 
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3 Deformations 



In this section, we consider infinitesimal deformations of a hypersurface S, and how the various geometrical 
quantities that characterize it change under deformation. A one-parameter deformation of the hypersurface 
is described by the functions X(£ a , u). An infinitesimal change of the embedding functions 

X -» X(0 + 5X(0 , (25) 

is characterized by the infinitesimal vector 5X(£) = <9„X(£ a , u)\ u= qSu. This vector can be decomposed 
into its components normal and tangential to the hypersurface E: 

5X = £lX + <5|| X = $n + $ Q e a . (26) 

The components 4>, <& a have the dimensions of length. We will assume that they are much smaller than 
any characteristic length in the system, such as the curvature radii. 

We first determine how the basis {e a ,n} changes under deformation. The key observation is that d u 
and d a commute, so that the change in the tangent vector satisfies 

5e a = d a (SX). (27) 

We now decompose <5X into its tangential and normal parts according to (|2^) so that 

6e a = (V a $ h )e b - K ab $ b n + (V„*)n + <5>K ab g bc e c , (28) 



where we have used the Gauss- Weingarten equations (15) and (|16|). We thus have for the deformations 
induced by 5Xm and 5'K±, respectively, 

5||e„ = (V a $ b )e fc - K ab ^> b n , 5 ± e a = (V a $)n + <PK ab g bc e c . (29) 

For the deformation of the induced metric, it follows that 

5\\g ab = V„$6 + V 6 *a, S ±9ab = 2K ab <P . (30) 

The tangential deformation is just the the Lie derivative along the surface vector field $ a . The normal 
deformation summarizes the geometrical content of the extrinsic curvature as one half of Lie derivative of 
the intrinsic metric along the normal vector field. Note that for the inverse metric one has that 5±g ab — 
—2K ab Q>. It follows from these relations that the first order deformation of the infinitesimal area element 
is 

5|| dA = dA V a $ Q , S±dA — dA K $ . (31) 

This last expression encodes the geometrical content of the trace of the extrinsic curvature as the relative 
change of area per unit normal deformation. 
The variations of the normal vector, 

5\\n = K ab <S> a g bc e c , S ± n = -(V a $) 5 ab e fc , (32) 



follow readily from the defining relations (11) and (|12|), together with (|29|). 

For the remainder of this section we will focus on normal deformations of the hypersurface. The reason 
is that infinitesimal tangential deformations correspond to reparametrizations of the surface, and, as shown 
in the next section, they contribute only boundary terms to the variation of global quantities. 

It is possible to evaluate the first order normal variation of the intrinsic scalar curvature 1Z either 
intrinsically or extrinsically, the latter using the Gauss-Codazzi equation. Let us consider the former 
approach. Using the definition of the Riemann tensor given by (j^), we have that the variation of the 
Riemann tensor is (see e.g. |24j) 

STZ a bcd = V c (6T a db ) - V d (6T a cb ) , 
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which implies 

6K ab = V c (ST c ab ) - V b (ST c ca ) . 
For the scalar curvature we have then that its normal variation is 

6 ± K = (6 ± g ab )K ab + g ab (5 x 1l ab ) 

= -2K ab K ab <S> + V c ( 5 afc <5 ± r^) - v a (5 ± r c ca ) . 

Now, an arbitrary variation of the metric induces a variation 

$Kb = \g cd {^b5g a d + V a Sg bd - V d Sg ab ) 

in the Christoffel symbols, so that using the Codazzi-Mainardi equations, we have that the normal variation 
of the Christoffel symbol is 

S ± T c ab = A Q c V b $ + A b c V Q $ - A afc V c $ + (V a A b c )$ . (33) 

As a consequence of the Codazzi-Mainardi equations, the (apparently unsymmetric) last term on the right 
hand side is, in fact, symmetric under the interchange of tangent indices. We determine then that the 
scalar curvature varies according to 

5±ll = -2K ab K ab <5> + 2V a [{K ab - 5 afc A)V b $] . (34) 

Note that for the densitized scalar curvature we have 

8xV9~K = V9(- 2 GabK ab <S> + 2V a [{K ab - g ab K)V b <S>] . (35) 

It follows that the variation of the integrated scalar curvature is a pure divergence in two dimensions, as 
expected from the Gauss-Bonnet theorem (jlOj) . 

Let us consider now the first order normal deformation of the extrinsic geometry. We have that 

S±K ab = -(S±n)-d a e b -n-d a (Sxe b ) 

= (V c $)e c ■ d a e b -n-d a [(V b $)n + $A' bc e c ] 

= r^ fc V c $ - 9 a V b $ + K ac K c b $ 

= -V a Vfc$ + K ac K c b <P. (36) 

This is a remarkably simple expression. Using the contracted Gauss equation ( ^l|) we can cast this 
expression in the alternative form, 

S±K ab = -V a V b $ + {KK ab - U ab )^ . (37) 

For the trace of the extrinsic curvature this gives (recall that 6±g ab = — 2A afc $), 

5 ± K = -A$ + (7e- A 2 )$, (38) 

where A := g ab \7 a \7 b denotes the laplacian on E, and we have used the contracted Gauss equation (^2|). 

As a check of the consistency of the expressions we have derived, one can verify that they imply the 
vanishing of the first order variation of the contracted Gauss equation, i.e. 5±(1Z — K 2 + K ab K ab ) = 0. 

In the calculations of higher order variations, an essential ingredient is the commutator of the defor- 
mation operator and the covariant derivative on E. This commutator allows one to express deformations 
of the covariant derivative of any geometrical quantity of interest in terms of the covariant derivative of 
the variation of such quantity. For instance, when acting on an arbitrary second rank tensor A b in E, it 
is given by 

[s X) v a ]A b c = -(s x r d ab ) A d c + (s x r c ad ) V , (39) 
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where 5±r^ b is defined in (33). Also for this case, it is a useful check of the validity of these expressions, to 
verify that the first order variation of the Codazzi-Mainardi equation (po|), which involves such commutator, 
vanishes identically. 

A useful expression is the commutator of the deformation derivative with the Laplacian, acting on an 
arbitrary function /, 



[6±,A]f = -2K ab t>V a V b f - 2K ab V a fV b <S> + KV a fV a <S> - (V a K )(V a /)$ . 



(40) 



4 First variation 

In this section, we apply the formalism for deformations we have developed to the systematic determination 
of the Euler-Lagrange derivatives of the geometrical invariants which appear in the hamiltonian for the 
strict bilayer couple model 

F[X] = a dAK 2 + (3M + fiA + PV . (41) 



The first term is the bending energy, with a the bending rigidity, the constants fi, P, /3 are lagrange 
multipliers that enforce the constraints of constant area, constant volume and constant area difference, 
respectively. For the sake of simplicity, we have not included the non-local bending rigidity term, necessary 
in a realistic description of lipid vesicles |^5|, ^(| . 

We first comment on tangential deformations. For the surface area element we have that under a 
tangential deormation it transforms according to (|3l|). In addition, any surface scalar satisfies S^f(X) = 
& a d a f(X). Thus, any reparameterization invariant functional of the form 

F[X] = J dAf(X) , (42) 

deforms tangentially as 

S {l F[X} = J dAV a [$ a f(X)}. (43) 

Using Stokes theorem this becomes 

5 {l F[X] = J ds Va <P a f(X), (44) 

where the integral ranges over the boundary of E and rj a denotes the unit normal to the boundary into E. 
If this boundary is empty, the integral vanishes identically. It is worthwhile emphasizing that in the case 
of a surface with a boundary, this integral in general will be non-vanishing, and it is no longer correct to 
neglect the tangential variations of the surface. In fact, in the variational principle, the vanishing of such 
terms will determine the boundary conditions to be imposed on the shape functions (see e.g. p7|). 

If the surface has no boundary, the tangential part of the variation can always be associated with a 
reparametrization. Since we are interested, to begin with, in quantities that are invariant under surface 
reparametrizations, we disregard this contribution and focus on normal deformations of the geometry of 
E. We write 

S ± F[X] = JdA [£[f] $ + V Q l/ ( a 1} [/]} . (45) 

The divergence comes about when we integrate by parts to remove all derivatives from the normal defor- 
mation <E>. The subscript in the vector appearing in the second term refers to the order of the variation. 

We begin with the first order normal variation of A and V. For the normal variation of the total area 
it follows immediately from ( |3l| ) that 



S X A = dAK<5>. (46) 
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The difference in area between two surfaces separated by an infinitesimal constant normal distance $ is 
proportional to the integrated mean curvature. The area of the surface is extremal for arbitrary normal 
deformations, 8±A = 0, when K = at each point on the surface. 

One way to determine the normal deformation of the volume enclosed by the surface is to exploit the 
definition given in (13). We have 



S - V = N 



/ d N Z pxV?)n-X + V9(5xn)-X + Vpn-5 ± X] 
i-j- J dA[$Kn-X- (V *) 9 ab e 6 • X + $] . (47) 



N 

Integrating by parts the second term and dropping a total divergence gives 
5±V 



N 
N 



^_ J dA [$ K n • X + $ g ab d a e b ■ X + $ g ab e b ■ e a + $] 
— / dA [$ K n • X - $ K n • X + iV$ + $1 

$ , (48) 

where we have used the Gauss- Weingarten equations ([l5|) to obtain the second term of the second line. 
This final expression should come as no surprise. The infinitesimal change in volume is simply proportional 
to the area of the surface times the normal displacement. 

We now examine the integrated powers of the mean extrinsic curvature. The deformation of the density 
y/gK , using @, f|), is 

S ±y /g X = -A$ + $ . (49) 
We thus find for the total mean curvature, up to a total divergence, 

S±M = I dAK<P. (50) 



It is interesting that this expression depends only on the intrinsic geometry of the surface. If the scalar 
curvature vanishes, so also does S±M . This is to be expected in a two-dimensional surface where 1Z = 
implies that it is flat. It is, however, a non-trivial statement for higher dimensions. 
For the variation of the integrated second power of the extrinsic curvature we have 

S ± JdAK 2 = J dAK [-2A$ + (211 - K 2 )<P] ; (51) 

We now integrate by parts twice to obtain, again up to a total divergence, 

6± JdAK 2 = J dA [-2AK + (211 - K 2 )K] $ . (52) 

We are now in the position to derive the Euler-Lagrange derivative for the model defined by the Helfrich 
hamiltonian @. Using the expressions ©, @, ©, ©, in @, we find 

£ = —2aAK + (2aK + /3)1Z — aK 3 + XK + P . (53) 

At equilibrium, £ = 0, it is known as the shape equation. Note that it involves four derivatives of the shape 
functions. It determines the equilibrium configurations of lipid membranes described by the hamiltonian 
(|l|). Recently it was shown how to cast the shape equation in the form of a conservation law [E8| . 
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5 Second variation 



In this section, we derive the second variation of the hamiltonian (^Tj). Before we proceed, we note that 
because the shape functions X are the variables to be varied, we have that S±& ^ , to second order, and 
similarly for the remaining three possibilities, 8±Q a , 6\\$> and 8\\& a . In fact, using (p9j), (j32|), we have 

g ± $= -$ a ■ v a $ , 8\\$ = K ab <S> a <i> b , (54) 

together with 

6 ± $ a = 3>V a <P + K a b <i> b <S> , 8^ a = $ fc V Q $ b -K a b <P b <S>. (55) 

This appears to suggest that the computation of the second variation is a formidable task and that the 
decompostion of deformations into tangential and normal parts will not be so useful. As we will show, 
the formal apparatus we have developed simplifies the calculation enormously and all of the troublesome 
tangential terms are collected in a divergence when the Euler-Lagrange equations are satisfied: no actual 
error is incurred in setting 8±$ = and neglecting $ a in the sequel. 



Let T = Jgf, where / is a scalar, introduced in (42). At first order, we have 

SF = Jd N ^ [6\\F+5 ± r] = J d n £ yg£® + V a {V {1) [T]+T$ a )] , (56) 

where VftJJF] = y/gVft)[f] i s the densitized vector defined by Eq.([45|). Now the divergence of a vector 
density is independent of the surface afirne connection r^ fc , and thus variation of the divergence of a vector 
density V a is equal to the divergence of the variation, 

5(X7 a V a )=X7 a (5V a ). (57) 

We thus have at second order, 

8 2 F = Jd N £S {^£<P} , (58) 

modulo another divergence. Thus, whereas for the second variation of T we cannot legitimately discard 
the divergence, we can for F. Furthermore, we have 

S Wg£$] = 8 ± (^/g£)$ + V a (^g£t> a )<P + ^g£8<S> . (59) 

Thus, modulo the Euler-Lagrange equation, 8 [y/g£$] = 8± (yfg£)Q. But 8± (^Jg£) can always be expressed 
as £<&, for some local differential operator £, and therefore the second order variation of the hamiltonian 
can itself be written in the form 

8 2 F[X] = [ dA<P £<$>. (60) 



It must be remembered, however, that this is not true for the individual terms contributing to F, for which 
£ ^ 0. When the second variation of these terms is written down it, it should be understand that the 
expression refers to the second normal variation; off diagonal terms ~ as well as purely tangential 
terms ~ $ a &b have been suppressed. 

A second order variation can also be seen as a different deformation — > -X" M (£) + <I>'n applied to 
the first order variations, and then letting $' = $. From a computational point of view, this is clearly 
equivalent to a repeated application of the normal deformation operator 8± . 

For the second order variation of the enclosed volume no work is needed, since, using (|3l|), it is simply 
related to the first variation of the area with 

8 2 V = 8 / dA$= / dAK<$> 2 . (61) 
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In order to derive the second variation of the area, various strategies are possible. On one hand one 
can simply calculate the second variation of the metric. Using ( |37j ) one obtains, 

Slg ab = 2 {S±K ab ) $ = -2$ [V a V b $ - K a c K cb ] $ , (62) 

so for the area element this gives 

siVg = VH-^ + T^ 2 ), (63) 

and the second variation of the area is given by the well-known expression (see e.g. 

S 2 A = J dA$ {-A + K} $ . (64) 

There is an alternative route to the second variation of the area, which becomes increasingly useful as 
one considers higher variations. The idea is to exploit the natural hierarchy in the variations illustrated 
by the relations 

Using (pf9|), this gives directly (|63|) . The usefulness of this approach is apparent: the n th variation of the 
surface volume element can be expressed in terms of the (n— l) th variation of the densitized mean extrinsic 
curvature. 

We now exploit this observation to derive the second variation of the densitized mean extrinsic curva- 
ture. We rewrite ( p0| ) as 

MVPO = Vs^ $ + VaVf 1) [if], (65) 

where we have defined the vector density VSJ-fT] := —y/gV a §. Using the fact that the variation of the 
divergence of a vector density is equal to the divergence of the variation, we have 

= 2^g<S> [(K ab - Kg ab )V a V b <S> - G ab K ab $] + V a V ( Q 2) [K]) , (66) 
where we have used (|35|), and defined the vector density 

V ( a 2) [K] = 5±Vfa [K] = ^(2K ab - g ab K)\7 b <P . (67) 
We have therefore for the second variation of M, 

S 2 M = JdA {2$ [(K ab - Kg ab )V a V b - Q a bK ab \ $} . (68) 

Note that if the surface is flat then the second order variation, like the first, vanishes. Thus, in a gaussian 
approximation about a flat background geometry, this term is absent. 

Our final task in this section is to compute the second variation of the densitized squared mean extrinsic 
curvature. Unfortunately, this quantity does not fit naturally in the hierarchy of variations we have been 
exploiting. One approach is to calculate the second variation of K directly. From (^8|), using (HJ), (ffO[), 
we have that 

5\K = 4X afc $V a V 6 $ + {2K ab - A; 9 ah )V a $V b $ + 2K(K 2 - TZ)^ 2 - 2TZ ab K ab ^ 2 . (69) 

We prefer to capitalize on the relative simplicity of the expressions for the variation of the densitized mean 
curvature and use the identity 

SKVdK 2 ) = 2K Sl(VgK) + 2[MVW] V 1/2 - 5K 2 <PS ± (^g-K) + 2^K^ 2 . (70) 
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A short calculation gives 

8\(yfgK 2 ) = ^[2(A$) 2 + A<S>KK ab V a V b <5> + (K 2 - 4^)$A$ 

+ (2ft 2 + 2K A - 5KK 2 - Ag ab K ab K)<P 2 ] + 2if V a V ( a 2) [K] , (71) 

where V" 2 ) [X] is defined by ( |S7|) . We isolate a total divergence in the final term using 

2XV a V ( Q 2) [K] = -2(V a X)V ( a 2) [X] + V a V ( a 2) [K 2 ] , (72) 

where we define 

V ( Q 2) [X 2 ] = 2y/g [A$V a $ - $V a A$ + $(2K ab - Kg ab )Vb<f>] , (73) 

to obtain 

S 2 J dAK 2 = J dA $[2A 2 + 4KK ab V a V b + (K 2 - 4K)A + 2K 2 

- 2(V a K)(2K ab - g ab K)V b + 2K i -5KK 2 -4g ab K ab K}<P, (74) 
We can simplify this expression using 

-2 J dA${\7 a K){2K ab - g ab K)\7 b $> = J dA[2K ab \7 a S7 b K - KAK + (V a iT)(V a X)]$ 2 , (75) 

to arrive at the final expression 

S 2 J dAK 2 = J dA<5> [2A 2 + 4:KK ab V a V b + (K 2 - 411) A + 2K ab V a V b K 

- KAK + ( V a K) ( V a K) + 2V? + 2K 4 - bKK 2 - 4g ab K ab K] $ . (76) 



Therefore, for the second order variation of the Helfrich hamiltonian ( |4l[ ), as expressed in (J60|), we can 
write the local differential operator C as 



C = 2aA 2 + 2AA + 2A ab V a V b + 2B , (77) 



where 



2A = a(K 2 — 41Z) - 2(3K — A , (78) 

2A ab = 4aKK ab + 2[3K ab (79) 

2B = al2K 2 + 2K 4 -5KK 2 -4g ab K ab K + 2K ab \7 a V b K 

- KAK + (V a K) ( V a K)} ~ 2f3g ab K ab + XJZ + PK . (80) 

We note that A ab is symmetric. Furthermore, the operator C is always self-adjoint, by which we mean 
that 



dA$i£$ 2 = / dA<$> 2 C<5>i ■ (81) 

Thus, the eigenvalues of C are assured to be real valued. To see this, consider the terms that could 
potentially spoil self-adjointness originating in the contributions proportional to a with two derivatives of 
$. Let A ab = A ab + g ab A. Now, 

$i afc V a V 6 $ = $V Q (i afc Vf,$) - $V a i ab V 6 $ 

= $V Q (i afc V 6 $) + iv a V 6 i ab $ 2 , (82) 
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modulo a total divergence which we discard. In this form the operator A ab V a V& is manifestly self-adjoint: 
we have 

*iV (A B6 Vi$ 2 ) - $ 2 V 6 (i ab V a $i) 
modulo two total divergences. We can write 

C = 2aA 2 + 2W a A ab W b + 2B 1 , (83) 

where 

2B X = a[2K 2 + 2K 4 - 5KK 2 - \Q ab K ah K + 2A(K 2 - 11)} - 2f3Q ab K ah + X1Z + PK. (84) 



We note that the expression Eq.(76) agrees with the one obtained in |16|, |17|] , in the terms containing 
derivatives of $. These previous studies were focused on the effect of short wavelength fluctuations, 
where only these terms contribute. We also note that the expression we have obtained disagrees with 
the one obtained previously in fl2|| . Agreement is restored in the special case of spherical and cylindrical 
configurations. 

It is clear that a rigid translation of the surface should not alter the second variation, at equilibrium. 
This provides a non-trivial check of our expression (77). As it is not entirely obvious, we will outline the 
details. Throughout, we will use extensively the Gauss- Weingarten equations (|l^), (|l6|), and the contracted 
Gauss-Codazzi-Mainardi equations (pit), (p2]), (23). Consider a constant normal deformation (5X = a, such 
that a • e = 0. We have for the derivatives of $ that appear in (|77j), 

V a V fc $ = -K ac K c b (n-a), 

A 2 $ = [-2K ab V a ^ b ^~^ a K){V a K)-A{K ab K ab ) + (K ab K ab ) 2 ]{n-a). 

If we now insert these expressions in the second variation, we find 

S 2 F = J dA {a[-2K ab V a S7 b <S> - (V a K)(V a K) - A(K ab K ab ) - KAK + 2TZK 2 - K 4 } 

+ (3KK + \K 2 + PK}{n-a) 2 . (85) 

We now integrate by parts the terms that involve derivatives of the extrinsic curvature. The desired 
simplifications follow from the crucial identity, for any surface vector V a , 



J dA (V a y a )(n-a) 2 -0, 



up to boundary terms. This is a consequence of V a (n • a) = K ab (e a ■ a) = 0, by hypothesis. Using this 
identity, we easily arrive at 

S 2 F = J dA K£(n ■ a) 2 , (86) 

where £ is the Euler-lagrange derivative defined by (|53|). Therefore at equilibrium the second variation 
vanishes for rigid translations. 

Note that a similar check has been used for the second variation of particular configurations, as, for 
example, spheres in [O, M|. However, not all of the terms appearing in the second variation can be checked 
considering only symmetric configurations, since many terms simply vanish identically in these limits. 



6 Concluding remarks 

In this paper we have presented a covariant geometric approach for examining the variation of geometric 
models of lipid membranes. For concreteness, we have restricted our attention to the first and second 



variation of the rigid bilayer couple hamiltonian (41). We note that it is straightforward to specialize the 
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second variation in the form (|77| ) to axisymmetric configurations. Moreover, our approach can be extended 
straightforwardly to other geometric models. However, computational difficulties are to be expected if 
the membrane has a free edge and it is no longer legitimate to throw away boundary terms. When 
considering higher order variations, a different difficulty arises. It is no longer justified to neglect tangential 
deformations. We will address this issue in a future publication. 
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